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We characterize all maximally entangling bipartite unitary operators, acting on systems A, B of 
arbitrary finite dimensions cLa < ds , when use of ancillary systems by both parties is allowed. Several 
useful and interesting consequences of this characterization are discussed, including an understanding 
of why the entangling and disentangling capacities of a given (maximally entangling) unitary can 
differ and a proof that these capacities must be equal when oIa — ds- 



I. INTRODUCTION 



A key question in quantum information theory is to understand the communication capabilities of quantum 
channels, wherein information is encoded in a quantum system which is then sent through the channel, which 
will generally introduce noise into the state of the transmitted system. A noisy quantum channel can be 
modeled as a unitary interaction between the system and its environment, and it is common to assume 
*£h that the environment starts out in a fixed pure state. One can, however, imagine a more general scenario, 

where the system and environment are allowed to have a completely arbitrary initial state. Then, we are 
considering the action of a unitary gate between two systems, and by varying their initial state, we can 
seek to maximize the amount of information communicated between the corresponding parties, where this 
communication may be in the form of classical information, quantum information, or both. In the case of 
quantum information, there is a close relationship to the amount of entanglement that can be produced 
by the given unitary interaction, and this is the question of interest to us here: What is the capacity of a 
bipartite unitary gate to generate entanglement [l|? 

We consider a unitary U acting on systems A and B held by Alice and Bob, respectively, system A 
pT) ' described by Hilbert space Ha and B by Hb, these Hilbert spaces having dimensions 6a < ds- Alice and 

f-*- , Bob are allowed the use of ancillary systems, a held by Alice (H a ) and b held by Bob (Hb)- It is well known 

that the use of ancillary systems increases the capacity of a unitary to generate entanglement 0, Q ■ If the 
input state on AaBb is \&i n ) and the output is then l^out) = Ia®Ib®M\^in), with I a (b) the identity operator 
00 ' on a (b), the capacity to generate entanglement is denned as 



£(U) = sup (E(^ out ) - E(% n )) , (1) 

|*in> 



where E(^) measures the entanglement of |^). The maximum possible value of £ (U) is 2\ogd,A, since any 
U can be simulated by LOCC using this amount of entanglement as a resource (the state of Alice's system 
can then be teleported to Bob and back) and LOCC cannot increase entanglement [4{. In this paper, we will 
only be interested in those U that are maximally entangling, that is, those that can increase entanglement 
by 21ogd J 4 ebits with some choice of \^i n )- 

In general, it is not known how large the ancilla need be to maximize the generation of entanglement for 
a given U, a significant barrier to understanding the entangling capacity of unitary interactions. However, 
in the case that U is maximally entangling, it has been shown that one can restrict consideration to d a = dA 
and db = ds [![■ There it was also shown that in this case, one may use an initial state that is product, 

\V in) AaBb = \$>} Aa ® |*)b6, (2) 

with \&/Aa = Sfc=i \fy a\k) a I Vd\ a maximally entangled state. In the next section, we use these results 
to characterize all maximally entangling unitaries for any dimensions dA,ds- Then, in section IIII1 we 
deduce several consequences of this characterization. Finally, in section IIV1 we summarize what has been 
accomplished. 
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II. CHARACTERIZATION OF MAXIMALLY ENTANGLING UNITARIES 



Our goal is to establish a characterization of maximally entangling bipartite unitaries. To that end, we 
will find it convenient to expand U, assumed to be unitary and maximally entangling, in terms of a finite 
group G, elements f,g€ G, group multiplication represented by fg. Thus, we have 

u = j2nf)^w(f), (3) 

fee 

where W(f) act on H B , {r(/)} are a set of unitary matrices forming a representation of G and acting on 
Ha, r(/)r(g) = g)T(fg), with the quantities g) constituting a factor system for which we have 
that \fJ.(f,g)\ — 1 V/,<7 due to the fact that T(f) is unitary for each /. When fJ,(f,g) — 1 V/,<?, we have an 
ordinary representation of the group; otherwise it is known as a projective representation. 

We know that such an expansion is always possible, since there exist groups of order \G\ = d\ that have 
representations forming a complete basis of the space of cLa X 4a matrices; the generalized Pauli operators 
provide one such example. However, for many unitaries, smaller groups are certainly possible, so we need 
to address the question of how to choose G. In Q, we used this type of expansion of bipartite unitaries to 
develop protocols for implementing 11 using local operations and classical communication (LOCC) with prior 
shared entanglement as a resource. For a given IA and any G with which such an expansion of U is possible, 
we showed how to deterministically simulate U by LOCC with a resource state having entanglement equal 
to log \G\. Since LOCC cannot increase the entanglement, we see that log \G\ must be at least as large as 
the amount of entanglement that hi can generate. This means that for the maximally entangling unitaries 
we are considering here, which have the ability to generate 2 log cU ebits, we need a group of order \G\ > d\. 

Writing the (to this point unknown) initial state on Bb as 

d B 

I*)s6 = E M mn\n) b \m) B , (4) 

m,n— 1 

the action of hi on the input state \^i n ) of @ yields 

\^out) = ^[i" a ®r(/)]|$) Aa ® [h®W{f)]\^) Bh 

/6G 

d A 



\= E I^U^E™^ E M mn \n) b W(f)\m) B 



^ A j,k=l feG m,n=l 

dA 

= J- E \ k )o\j)A ® \b jk ), (5) 

where [r(/)] jfc is the jk matrix element of T(/), this basis chosen for convenience to be that which completely 
reduces the matrices into irreducible representations (the finest block-diagonal form of these matrices). 
We have defined 



IV) = V^E[ r (/)U E M mn \n) b W(f)\m) B . (6) 

/GG rn,n— 1 

Assuming that \&i n ) is an optimal input, achieving the maximal entanglement generation of 21ogcU ebits, 
we see immediately from ([5]) that the states \bjk) must form an orthonormal set, Sjj'Skk' — (bj'k'lbjk)- This 
implies, first of all, that for each fixed j, k, 3f such that [r(/)] J 7 C ^ 0. Recalling that we have chosen the j, k 
basis to be that which completely decomposes matrices T(/) into irreducible representations, we see that 
these matrices are themselves an irreducible representation for G of dimension cLa- Therefore, the choice 
of G is restricted to one which has an irreducible representation of this dimension. Since the sum of the 
squared dimensions of all irreducible representations of G is equal to |G|, we here have another (related) way 
of seeing that \G\ > d\. As mentioned previously, we can always choose a representation by the generalized 
Pauli matrices, for which \G\ = d\ (with this choice we have a projective irreducible representation, and 
for the given factor system, this is the only irreducible representation for G), and we will assume this 
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choice has been made throughout the remainder of this paper. It is often convenient to define the 
such that the factor system satisfies /i(e,g) = fi(g,e) = fi(g,g 1 ) = 1, V<? G G. One possibility is to use 
r(/) = T(m, n) = e ie ™X m Z n , with 9 mn = n[mn (mod d A )]/d A - 

It is shown in appendix|X]that as a consequence of Schur's orthogonality relations for group representations 
, the orthonormality condition on states \bjk) is equivalent to a corresponding orthonormality condition 
on operators W(f), 

Tr [W(f)MM^W( g y] = ±5(f, g), V/, g € G, (7) 

where S(f,g) — 1 when f — g, and otherwise is equal to zero. Thus we have our main result: 

Theorem 1. The bipartite unitary IA is maximally entangling iff there exists a positive semi-definite 'metric' 
MM^ such that ([?[) is satisfied V/, g € G, where operators W(f) are obtained from an expansion of IA as 
in (0), with the taken to be the generalized Pauli operators. The operator M defines an optimal input 
state on systems bB through |^). 

In the next section, we discuss consequences of this result. 



III. CONSEQUENCES OF THEOREM [TJ 
Consequence 1. Method to check if Li is maximally entangling. 

Given bipartite unitary IA, theorem [T] provides a method of determining whether or not IA is maximally 
entangling. One need only expand IA in terms of the generalized Pauli operators, identify the set of operators 
{W(f)}, and then check to see if there exists a positive semi-definite operator to play the role of MM^ such 
that ([7]) is satisfied. One way to do this is to form all products, W(g)' ! W(f), f ^ g, reshape each into a 
column vector (such as by stacking individual columns of each product one on top of the other) and collect all 
these columns into a matrix. The nullspace of this matrix corresponds (by reshaping vectors in this nullspace 
back into matrices) to the space of all operators orthogonal to the W(gy\V(f), f ^ g, as is required to 
satisfy ([7]). One then needs to search, perhaps numerically, for positive operators in this nullspace. This is 
relatively easy to do, at least for small enough nullspaces. 

Consequence 2. Design of maximally entangling unitaries. 

Theorem [1] also allows one to design unitaries that are maximally entangling. This amounts to choosing 
operator M and a set of d\ linearly independent operators W(f) that satisfy ([7]). In addition, there is also 
the necessity that the chosen set of W(f) are such that IA is unitary. When the dimensions are not too large, 
it is straightforward and reasonably fast to numerically generate a maximally entangling unitary in this way 
(for d A = 4, rig = 8 it takes less than 10 minutes on my laptop). 

Consequence 3. Characterizing maximally entangling interaction Hamiltonians for two-qubit systems. 

A characterization of two-qubit maximally entangling Hamiltonians H has been given in 0. Using the 
well-known result [2[ that up to local unitaries, every two-qubit unitary may be written as IA = e~ lH (the 
usual factor t/H is here absorbed into the definition of H for notational convenience) with 

H = a i a i® a h ( 8 ) 

j=x.,y.z 

they showed that for IA to be maximally entangling, it must be that cos 2 a x = 1/2 = cos 2 a y , with the value 
of a z being unconstrained (permutations of {x, y, z} are also allowed, of course). In appendix iBl we provide 
an alternative proof of this result based on ([TJ. This means there is a continuum of maximally entangling 
two-qubit unitaries ranging from the double CNOT (a z — 0) to the SWAP (cos 2 a z = 1/2). 

Consequence 4. Operators W(f) must form a linearly independent set. 
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This is easily proven, as is shown at the end of appendix [XJ Notice also that by theorem 4 of Q and for 
whatever group G and representation T are chosen for the expansion of maximally entangling IA, the number 
of linearly independent operators in the collection {r(/)} is d\, because only the single 44-dimensional 
irreducible representation appears in these matrices. This is consistent with the fact that the Schmidt rank 
of IA must be at least as large as the ratio of the Schmidt rank of the output state to that of the input state. 
That is, since our input state has Schmidt rank of one and the output state has Schmidt rank of d\, IA must 
have Schmidt rank of d\ as well. 

Consequence 5. Input state on Bb is uniquely determined up to local unitaries when dA — du, and must 
be a maximally entangled state. 

This was proven in Q ; we provide an alternative proof based on ([7]) in appendix [X] 
Consequence 6. Why the entangling and disentangling powers can be unequal. 

It is now easy to see for a maximally entangling unitary how the entangling and disentangling powers can be 
unequal Recall that the disentangling power of IA is just the entangling power of W . Therefore for the 
disentangling power, we must replace the set {W(/)} by in ([7]). Then, for IA to be maximally disen- 

tangling, we require the existence of an M'M'^ orthogonal to the set of operators {W{g)W(f)^ }, V/ 7^ g € G, 
whereas for maximally entangling, the orthogonality requirement applies to the generally different set, 
V/ ^ j £ G. In addition, there is the normalization condition for f = g, and this again 
applies to a generally different set of operators in the two cases. As an example, Q provided the origi- 
nal demonstration that the entangling and disentangling powers can be unequal by constructing a specific 
maximally entangling IA and then showing that lA^ has strictly less than the maximum entangling power. 
We have calculated the W(f) for their IA and find that it is easy to satisfy (J7J with these W(f) (set 
MAft = + |3) B (3|]/2), but find (numerically) that it is not possible to do so when the set {W(f)} 

is replaced by {W(/) t } (one choice that almost works is to set M'M't = c |l) B (l| + c(|2) B (2| + |3) B (3|), 
which satisfies orthogonality, but the normalizations cannot all be the same no matter how cq, c are chosen). 

Consequence 7. Entangling and disentangling powers are equal for maximally entangling unitaries on dxd 
systems. 

It was shown in Q that the entangling and disentangling powers of any IA are equal when dA — 2 = c?b. 
We can now extend this result to arbitrary dimensions dA = <^b when restricting to maximally entangling 
unitaries. From consequence [5l we have that MM' must be proportional to Ib- Therefore, a replacement 
{W(f)} — > {W(fy} makes no difference whatsoever in ([7]), from which this claim follows immediately. That 
is, when dA = ds and IA is maximally entangling, then IA^ is also maximally entangling. 

Consequence 8. If ds is large enough compared to dA, it can be that no ancillary system is needed on 
Bob's side. 

We here provide a construction of operators W(f) corresponding to IA for which system b is not needed. 
This requires only that the first columns of the different W(f) operators are mutually orthogonal and have 
norm equal to 1/ 'dA (the remaining part of each W{f) is unconstrained apart from the requirement that IA is 
unitary). Then we have that the matrix element (l\W(g)^W(f)\l) = 5(f,g)/d 2 A . Choosing MM* = 
shows that J7J) is satisfied V/, g. This choice of MM^ corresponds to a product state across B/b, so system b 
never plays a role and may be discarded. Recalling that there are d? A different W(f) operators, the mutual 
orthogonality of their first columns is possible only when the length ds of those columns is at least d\. 
Hence, this construction is only possible when ds > d\. Then there is a (^-dimensional subspace of 'Kb 
that becomes maximally entangled with systems Aa, the remaining space not being involved in the process. 
Thus, it is almost as if system B has the ancillary system already embedded within itself, which is most 
clearly understood when ds = d\ = dA x dA- In this case, B can be thought of as itself consisting of two 
<iyt-dimensional systems, one of which plays the role of ancillary b. 

IV. CONCLUSIONS 

We have given a characterization of all maximally entangling bipartite unitaries for any dimensions dA < 
ds- This allows one to check if a given unitary is maximally entangling, to construct maximally entangling 
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unitaries, and to determine optimal input states that achieve the maximal generation of entanglement. It 
also provides an understanding of why the entangling and disentangling capacities can differ, as well as a 
proof that this can only happen when ds > d A . We also saw that for ds > d\, it is possible that no 
ancillary system is needed on Bob's side. Finally, we have given an alternative method of characterizing 
maximally entangling Hamiltonians for two-qubit systems Q. An interesting open question is to determine 
what Hamiltonians can be maximally entangling in higher-dimensional systems. 
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Appendix A: Proof of theorem [T] 

Here we show that orthonormality of the states \bjk) defined in |6]) is equivalent to condition ([7]) on 
operators W(f), which appear in an expansion of U of the form ([3]) with \G\ — d\. From |6j), we have 

(b fk ,\bik) = d A ]T [T(g)]*, k , [T(f)] jk £ £ M^ n M mn {n'\n){m'\W{g)^W{f)\m) 

f,g£G m,n— 1 m' ,n' — 1 

= d A [r(s)]*^r(/)] ifc J2 [MM^ mml {m'\W( g yW(f)\m) 

f,g£G m',m— 1 

= d A FWv Tr [MM^W( g yW(f)} . (Al) 

f,geG 

First notice that if Tv [W(f)MM^W(g) i '] — 8(f,g)/d A , the right-hand side of this equation becomes 
Ylf \P(f)]j'k' [r(/)]jfc /4a- However, considering the d\ vectors jj k , j, k = l,...,d A , whose components 
(labeled by / £ G) are given by 



tf jk ) f = -j= [T(f)] jk , (A2) 



1 

fa 

then by Schur's orthogonality relations for irreducible representations Q and the fact that the rep- 
resentation is irreducible, these vectors form a complete orthonormal basis for the c^-dimensional space in 
which they lie (recall that \G\ = d\ is the dimension of these vectors). That is, 

E FifTj'k' i T (f)hk l d A = 6 jf 5 kk ,, (A3) 
/ 

which yields one of the implications we sought to prove. 

To prove the converse, define d\ x d\ matrix O, with matrix elements labeled by /, g € G given by 

[0] gf = Tr [W(f)MM*W( 9 y] . (A4) 
Then if (bj> k >\bjk) = SjfS kk >, (IA1|) can be written as 

— S jf 8 kk , = jt, k , ■ O ■ j jk . (A5) 

By (|A5|) . O ■ jj k is orthogonal to every vector in the complete basis of the 7- vectors except for one, that 
being jj k . Therefore, Vj, k, O ■ jj k is proportional to jj k , and the proportionality constant is equal to l/d^, 
independent of j,k, again by (| A5|) . Thus, we have that O — I/d\, where / is the d\ x d\ identity matrix. 
Finally, recalling the definition of O in (|A4j) , we have 

±5(f, g) = Tr [W(f)MM^W(g^] , (A6) 
a A 



G 



which completes the proof. ■ 
A necessary condition for (|A6|) to be satisfied is that the collection of \G\ = d\ operators W(f) are linearly 
independent. This is easily seen by contradiction, so assume they are linearly dependent. Then, 

0=$>(/W/), (A7) 

fee 

for some coefficients c(f) not all equal to 0. Multiply this expression by MM^W(g)^ for each fixed g € G 
and then take the trace to obtain from (IA6I) that 



= £ c(/)Tr [W{f)MM^W{gf] = (A8) 
fee A 

assuming (|A6I) . This says that c(g) = Wg G G, which contradicts the assumption of linear dependence and 
proves the claim. 

We now give an alternate proof (see also |5J ) that p = MM^ is uniquely determined when U is maximally 
entangling and cIa — ds- Indeed, by contradiction, assume both p and p' serve our purpose. Then from 
(TA61) . 

= Tr [W(f)( P - p')W(g)i] V/, geG, (A9) 

which must hold even when / = g. This says that for each f,g g G, W(f)(p — p') is orthogonal to W(g). 
However, as we have just seen, the d\ = d 2 B operators W(g) are linearly independent, hence span the entire 
space B(J-Lb) of operators acting on T~Lb- Therefore, it must be that 

W(f)(p-p') = (A10) 

for every f £ G. Now, choose coefficients e(/) such that 1b = 2/ e (/)^(/)> which can always be done 
since W(f) are a basis of B(Hb)- Multiplying (|A10[) by e(/) and summing over / we obtain = p — p', 
proving the claim. 



Appendix B: Two-qubit maximally entangling Hamiltonians 

Using (|8|) gives U = e~ lH — J2f kf<Jf ® aj with f = e,x,y,z labeling the group element (e is the identity 
element). From this we identify Wf — kf<Jf (cr e = /, the two-by-two identity matrix) , where 







SxSy^z i 


k x 


CxSyS z 


Sx C z j 


ky 


— S x CyS Z 


&x Sy Cz ■> 


k z 







and we've used the abbreviations Cf = cos a/ and sr = sin a/. Applying the condition (J7]) with MM' =1/2 
(because c?^ = rig), the orthogonality conditions (/ 7^ g) are automatically satisfied because the Pauli 
operators are themselves mutually orthogonal. Therefore, we only need to worry about normalizations 
(/ = 5 in (0), which give 



222, 

C x C y C z ~l 


2 2 2 
S x S y S z 


-1/4, 


2 2 2 
c x s y s z 


Us 2 2 c> 2 
r s x c y c z 


= 1/4, 


111 
s x c y s z 


,222 
r c x s y c z 


= 1/4, 


2 2 2 
S x S y C z 


,222 
r C x C y S z 


= 1/4. 



It not too difficult to show that these lead to the necessary and sufficient condition that two of the a's must 
have squared cosines equal to 1/2, the third a being unconstrained, which is what we set out to prove. 
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